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CLAIMS 

The embodiments of the invention in which an exclusive property or privilege is claimed are 
defined as follows; 

1-20 (previously cancelled) 
21-29 (cancelled) 

30. (New) A method of forming/defining and solving a model ofjk€ power network to eSfeff 
control of voltages and power flows in a power system, comprising the steps of: 
obtaining on-line/simulated data of open/close status of switches and circuit breakers in 
the power network, and reading data of operating limits of JfaT network 
components including PV-node, a generator-node where Real-Power-P and 
Voltage-Magnitude- V are given/assigned/specified/set, generators maximum and 
minimum reactive power generation capability limits and transformers tap position 
limits, 

obtaining on-line readings of given/assigned/specified/set real and reactive power at PQ- 
nodes, the load-nodes where Real-Power-P and Reactive-Power-Q are 
given/assigned/specified/set, real power and voltage magnitude at PV-nodes, 
voltage magnitude and angle at the reference/slack node, and transformer turns 
ratios, which are the controlled variables/parameters, 

initiating loadflow calculation with initial approximate/guess solution of the same voltage 
magnitude and angle as those of the (si ack/referen c^ node tor all the other nodes 
referred to as the slack-start, ***** °™ *A^*^ 

formivg and storing factorfeed gain matrices [Y0] and [YV] using the same indexing and 
addressing information for both as they are of the same dimension and sparsity 
structure, wherein said [Y0] relate vector of modified real power mismatches [RP] 
to angle corrections vector [A6] in equation [RP] = [Y8] [AG] referred to as P-0 
sub-problem, and said [W] relate vector of modified reactive power mismatches 
[RQJ to voltage magnitude corrections vector [AV] in equation [RQ] = [YV] [AV] 
referred to as Q-V sub-problem, 

restricting trans formation/rotarion angle <D P to maximum -48° in determining transformed 
real and reactive power mismatch as, 
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performing loadflow calculation by solving a Super Super Decoupled Loadflow mode l of 
the power network defined by set of equations [RP] = [Y9] \AQ) and [RQj = [YV] 
[AV] employing successive (16, IV) iteration scheme, wherein each iteration 
involves one calculation of [RPJ and [AO] to update voltage angle vector [0J and 
then one calculation of [RQ] and [AV] to update voltage magnitude vector [V], to 
calculate values of the voltage angle and the voltage magnitude at PQ-nodes, 
voltage angle and reactive power generation at PV-nodes, and turns ratio of tap- 
changing transformers in dependence on the set of said obtained-online readings, 
or givcn/scheduled/spccified/set values of controlled variables/parameters and 
physical limits of operation of the network components, 

evaluating loadflow calculation for any of the over loaded power network components and 
for under/over voltage at any of the network nodes, 

correcting one or more controlled parameters and repeating the calculating, performing, 
evaluating, and correcting steps until evaluating step finds no over loaded 
^components and no under/over voltages in^h^power network, and 
j?jfeet*Hg a change in the power flowing through network components and voltage 
magnitudes and angles at the nodes of J&(£ power network by actually 
implementing the finally obtained values of controlled variables/parameters after 
evaluating step finds a good power system or alternatively a power network 
without any overloaded components and under/over voltages, which finally 
obtained controlled variables/parameters however arc stored in case of simulation 
for acting upon fast in case the simulated event actually occurs. 
3 1 , (New) A method as defined in claim-2 1 wherein loadflow calculation involving formation 

and solution of super super decoupled loadflow model, employing simultaneous 
(IV, 10) iteration scheme is characterized in that it involve only one time 
calculation of real and reactive power mismatches in an iteration along with 
modified real power mismatch calculation, depending on super super decoupled 
loadflow model used, either by; 



m 



Rp p = [(APp/V„)-2:Gp q AVq]/ v p -for all nodes or (74) 
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Abstract - A generalized Super Super decoupled Loadflow 
referred to tw SSDL-X'X' method is presented along with its 
special case aionifc»t»tions leading to the proposed SSDI^YY 
method. Culn matrices of the SSDL-YY method are defined 
independent of rotation angles and they differ only in their 
diagonal terms. The diagonal terms in the gain matrix of Q-V 
suh-problem are characterized by the presence of transformed 
tipeeJftcd/k.nnwfi quantities. Rotation angles are limited to the 
maximum value of -4R\ and the slack-start as a guess solution is 
suggested for Increased efficiency. These along with simple 
modification of real power mismatches at PV-nodcs yields a 
robust and fast loadflow method. An advanced version of the 
Genera U purpose Fast Decoupled Loadflow (CFDL) U manifested 
as SSDL-BX method. The proposed method is tested and 
compared with Fast Super Decoupled Loadflow (FSDL), SSDL- 
UX* and CFDL methods under wide range of simulated system 
operating conditions and network parameters* The proposed 
method Is superior to the FSDL method under all most all 
conditions and can be more than twice faster than the SSDL-BX, 
and CFDL method* for systems having large R/X ratio branches. 

Indexing terms: Loadflow, Transformation, Decoupling 

T, NOMENCLATURE 

Y„ = Y M ^(cD w +180°) =G pq +jB pq : (p-q) th element of nodal 

admittance matrix without shunts 
y p ■ y p Z(0)^ ■ gp + jb p : total shunt admittance at any nodc-p 
Vp ■ Cp + jf P » VpZB p : complex voltage of any node-p 
V, ■ c» + jf p " V^6 B : complex slack-node voltage 
AOp, AV r 

+ JQr 



p 

AP 



n 

p - jAQ„ 
PSH. f jQSH r 



: voltage angle, magnitude corrections 
: net nodal injected power, calculated 
: nodal power residue or mismatch 
: modified nodal power residue or mismatch 
net nodal injected power, scheduled/specified 
C p " I ^Op"Cos<I>p+jSin<I>p:UnitaJ7rotation/u^nsfonnation 
m ; number of PQ-nodes 

k : number of P V-nodes 

n=m+k+l : total number of nodes 

q>p : node«q is connected to node-p excluding the case of q=p 

IT. INTRODUCTION 

Loadflow, in power system studies, is the most basic 
frequently performed steady state analysis of an electrical 
power network. Load flows are performed in system planning, 
operational planning* and operation control. Of the various 
methods proposed including the first time introduced in [1] 
decoupled loadflow referred to as Voltage Vector Method 
(WMX the Fast Decoupled Loadflow (FDL) [2] method was 
the siatc-of-thc-practice because of its simplicity, relatively 
high reliability, and computational efficiency. However, it is 
known to suffer poor convergence for systems having high 
R/X ratio branches. The General-purpose version of the Fast 
Decoupled Loadflow (GFDL) [3] is a simple modification of 
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the FDL with much improved convergence in the presence of 
high R/X ratio brancheSv However, it has been shown to have 
failed to converge to a solution when FDL converged [3]. 

The transformation or super decoupled approach (4), [5] is 
mathematically sound and therefore, more general and 
reliable. Rotation operators applied to the complex node 
injections and corresponding admittance values that relate 
them to the system state variables, transform the network 
equations such thai branch admittances appear to be almost 
entirely reactive. Better decoupling is thus realized. The initial 
methods of [4] and [5] were only slightly better convergent 
than those of [2] and [3], However, the super decoupling 
technique with gain matrices defined independently of rotation 
angles, restriction of rotation angle to maximum of -36° and a 
simple modification of rcaJ power mismatches at PV-nodcs led 
to the development of Fast Super Decoupled Loadflow 
(FSDL) [6] method. The process that began with the VVM 
leading to the state -of- the- practice FDL method and then to 
the current state-of-the-art FSDL method is being taken to the 
new state-of-the-art proposed Super Super Decoupled 
Loadflow (SSDL-YY) method in this paper. 

A. Decoupled Loadflow 

In a class of decoupled loadflow models, each decoupled 
loadflow model comprises a system of equations (1) and (2) 
differing in the definition of elements of [RP], [RQ], and [Y0] 
and [YV]. It is a system of equations for the separate 
calculation of voltage angle and voltage magnitude corrections. 

[RPJ = [Y0] [AO] (1) 

[RQ] - [YV] [AV] (2) 

8. Successive (IB, IV) Iteration Scheme 

In this scheme (1) and (2) are solved alternately with 
intermediate updating, Each iteration involves one calculation 
of [RP] and [A6] to update [0] and then one calculation of [RQ] 
and [AV] to update [V]. The sequence of relations (3) to (6) 
depicts the scheme. 



[AO] - [re] 1 im 

m = m + [AG] 

[AV^CYVT 1 [RQ] 
[VJ - [V] + [AV] 



(3) 

W 
(5) 

(6) 



The scheme involves solution of system of equations (1) 
and (2) in an iterative manner depicted in the sequence of 
relations (3) to (6). This scheme requires mismatch calculation 
for each half iteration because [RP] and [RQ] is calculated 
always using the most recent voltage angle and voltage 
magnitude values, and it is block Oauss-Seidel approach. The 
scheme is block successive, which imparts increased stability 
to the solution process, and it in turn improves convergence 
and increases the reliability of obtaining solution. 
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111. THE PROPOSED MODEL 

It con be seen from the Jacobian elements determined at the 
sJack-SHott for formulations of the power flow equations in tic 
nppend ix-A that Jacobian is symmetrical for 1=0, and 
unsym metrical for all other values of raised power 4 T\ Phase 
shifters ore assumed absent and shunt conductances arc 
ignored for the Jacobian to be symmetrical in the cose of 
formulation of ihc power flow equations in appcndix-A for 
1-0, for which the constant matrix model is given in [6]. 
Coupling between the active and reactive power flows in a 
power circuit is mainly attributed to the resistive component of 
the impedance. Decoupled loadflow models conscquendy have 
solution convergence problems in the presence of transmission 
lined ol r large R/X ratios. Better decoupling of the loadflow 
equations for networks having branch elements of high R/X is 
achieved through simple linear transformation of complex 
power injeelions as given in appendix-9.2 of [6], 

A, Super Super Decoupled Loadflow: Version: SSDL-X'X 9 

Transformed power flow equations (42) and (43) in 
append ix-B arc of the same form as the original real and 
reactive power flow equations contained in (32) and (33) in 
appendix- A, The only difference is that they involve 
conductances and susccptanccs values of the 
rotatcd/tro informed branch admittances. The most general 
Super Super Decoupled Loadflow model can be realized as 
SSDL-XTC, from which many versions arc derived. The 
elements of [RP], [RQ], [Y8J and [YV] constituting the model 
SSDL-X'X' arc defined by (7) to (19), 

RP r * [AP P ' + (GppVBpp-) AQp'yVp 2 :for PQ-nodes (7) 
RQp - [AQp f - (Gp,7e„,') AP p »]/% :fot PQ-nodcs (8) 
RP P * [APp/tKpV, 2 )] :for PV-nodes (9) 

Y9 N --I/X pq k and YV W = -1/X^' (10) 

and 



YV W = V + X-YV M 



GO 



Where, 

b^ - I [QSHp' - (CppV Bpp^PSHp'J/V, 2 - (2 - I) b p CosOp ( 1 2) 

APp 4 - APpCos<l> p + AQpSind>p :fer PQ-nodes (t 3) 

AQ P * = AQpCosOp -APpSincfrp :forPQ*nodes (14) 

PSry - PSH p CosO»p+ QSH p Sin<P p ;for PQ-nodcs (1 5) 

QSHp' * Q$HpCos<l> 0 - PSHpSina>p :for PQ-nodes (1 6) 

Cos<Dp - | [BppMcV +b pp 3 )] I 2 (°* to - 900 ) ( 1? ) 

Sinwp - - 1 (C^C^ *&pp*)] | a Sin (0* to -90*) (1 8) 

Kp- KBp^Bpp)! (19) 

Where. Gpp* and Bpp* arc the transformed values of Gpp and 
Bpp as given in the appendix-8 by (45) and used in (52) in the 
appcndix*C„ The factor K p of (19), determined experimentally; 
is restricted to the minimum value of 0.75, but its restriction is 
reduced to the minimum 0.6 if its average over all less than 1 ,0 
values at PV-nodes is less than 0.6. This factor is system and 
model independent, and does not require tuning. The factor *T* 
in (12) can take any value from I - - ©o 7 .... -2, -1,0, 1,2, ... f 



+ However* it has been determined experimentally that I = I 
gives the best possible convergence for all systems and it is 
also system and model independent. But in situations where 
the diagonal dominance in [YV] is influenced in power 
networks having capacitivc series branches and/or excessive 
shunt capacitivc compensation deteriorating convergence and 
reliability of obtaining solution, the best possible convergence 
can be obtained by tuning the factor "I" in (12). 

The definition of YBp,, in (10) is simplified because it docs 
not explicitly stale that it always takes the value of -B w for a 
brunch connected between two PV-nodcs or a PV-node and the 
slack-node. This fact should be understood implied in all the 
definitions of YOp, in this document. X^' is the transformed 
branch reactance defined in the appcndix-C by (51). 

#> Super Super Decoupled Loadflow: Version: SSDL-YY 

If unrestricted rotation is applied to complex branch 
admittance and transformed susceptonce is taken as admittance 
magnitude value with the same algebraic sign and transformed 
conductance is assumed zero in (44) of Appendix-B or (50) of 
Appcndix-C, the SSDL-XVX' model reduces to SSDL-YY. 
Though* this model is not very sensitive to the restriction 
applied to nodal rotation angles SSDL-YY presented here 
restricts rotation angles to the maximum of -48 n in (17) and 
(18) determined experimentally for Ihc best possible con- 
vergence from non-linearity considerations. However, it gives 
closely simitar performance over wide range of restriction 
applied to the nodal rotation angles say, from -36° to -90°. 



RPp = (APpCos0> p + AQpSintb p yv p * ;for PQ-nodcs 
RQ P = (AQpCosOp - APpSind>p )/Vp :fbr PQ-nodcs 



(20) 
(21) 



Y0 M - 



-Y 



:for branch r/x ratio £ 3-0 



- (B N + 0.9(Y PQ -B W )) :for branch r/x ratio > 3.0 

-Bpq :for branches connected between two PV- 
M nodes or a PV-node and the slack-node 



(22) 



YV = 



:for branch r/x ratio £ 3.0 
90VBP,,)) :fbr branch r/x ratio > 3.0 (23) 
Where, bp' of (1 2) becomes, 

bp I (QSH>Cos<P p - PSHpSin<Ppy V/ - (2 - 1) b p Cos<P p (24) 



Branch admittance magnitude in (22) and (23) is of the 
same algebraic sign as its susceptance. Elements of the two 
gain matrices differ in that diagonal elements of [YV] 
additionally contain the b p ' values given by (I I) and (24) and 
in respect of elements corresponding to branches connected 
between two PV-nodes or a PV-node and the slack-node. The 
model consists of (3) to (6), (20) to (24), (9), (1 1), (17), (18), 
and (19) with rotation angles restricted to the maximum of -4& A 
in (17) and (18). In case of systems of only PQ-nodes and 
without any PV-nodes, equations (22) and (23) simply be taken 
as YOpq - YV M ■ -Yp, , In two simple variations of the SSDL- 
YY model, one is to make YV m =Y8 pq and the other is to make 
YOp^YVp,. 

C. Super Super Decoupled Loadflow: Version: SSDL-XX 
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If no or zero rotation (<t> p »° 0*) is applied, the SSDL-X'X* 
model reduces to SSDL-XX, which is the simplest form of 
SSDL-X'X\ The SSDL-XX model comprises (3) to (6), (25) 
to (29). (9). (II), and (19). 

RP P - [AP P + (Opp/Bpp) AQ P ]/V P 2 rforPQ-nodes (25) 

RQp = [AQp- (Cpj/Bpp) APp]/V p :for PQ-nodcs (26) 

Y0 N ■ p 1 ,0/X,*, :for al I other branches 

«B M :for branches connected between two PV- 
L nodes or a PV-nodc and the slack-node (27) 

YV N - .1.0/X^ :for all branches (2S) 

Where. bp f of (12) becomes, 

b p ' - 1 {QSHp - (Gpp/Bpp) PSH P ]/V S 3 - (2 - T) bp (29) 

This is the simplest model with very good performance for 
distribution networks containing only PQ-nodes. The large 
value of the difference [(1/X M ) - B m ], particularly for high R/X 
ratio branches connected between two PV-nodes or a PV-node 
and the slack-node, creates modeling error in a symmetrical 
gain matrix model of a power network containing PV-nodes. 

/>, Super Super Decoupled Loadflow: Version: SSDL-BX 

if super decoupling is applied only to QV-sub problem, the 
SSDL-XX model reduces to SSDL-BX, which makes it 
perform better for systems containing PV-nodes» The SSDL- 
BX model comprises (3) to (6), (26), (30), (31), (1 1) and (29). 

RPp- APp/V n s :for all-nodes (30) 

Ye M =-3 M and YV pq = -l/X pq (31) 

It should be noted that Amerongen's General-purpose Fast 
Decoupled Loadflow model [3] has turned out to be an- 
approximation of this model. The approximation involved is 
only in (26) and (29). However, numerical performance is 
found to be better and more reliable than that of the 
Amerongen's method. This model can be referred to as 
advanced OFDL model. 

£ Stack-start 

The same voltage magnitude and angle as those of the 
slackmodc when used for all nodes as initial guess solution, 
provide Jacobian elements determined in the appcndix-A. With 
the specified magnitudes, PV-nodcs voltage magnitudes are 
adjusted to their known values after the first half iteration. This 
start procedure referred to as the slack-start, saves almost all 
effort of mismatch calculation in the first P-9 iteration as it 
requires only shunt flows to be calculated at each node. 

Ft Modd involving Global Corrections 

The non-linearity retaining technique [7] or an exact second 
order model [8] formulations for the solution of the loadflow 
problem in rectangular coordinates produce corrections to the 
initial guess solution in each iteration. However, the use of the 
corrections to the initial solution estimate in each iteration is 
reported first in [1] for the loadflow problem formulation in 
polar coordinates. A generalized technique for the formulation 



of all the loadflow models that produces corrections to initial 
guess solution in each iteration is given in appendix-D. It 
involves storage of the modified mismatch of the current 
iteration for adding into the modified mismatch calculated in 
the next iteration. 

G: Sequence of steps for the solution ofSSD/^YY model 

AM possible variations of the Super Super Decoupled 
Loadflow models arc derived from the most general model 
SSDLOCX' in the above and in [18]. It has been verified that 
SSDL-YY version is the simplest possible general purpose 
model that can perform well for power networks containing 
PV-nodes as well as those that contain only PQ-nodes. 
Therefore, detailed sequence of steps for the solution of only 
SSDL-YY model is given in the following, 

a. Read system data and assign an initial approximate solu- 
tion. If better solution estimate is not available, use the 
slack-start, 

b. Form nodal admittance matrix, and Initialize iteration 
count ITRP = 1TRQ= r = 0 

c. Compute Cosine and Sine of nodal rotation angles using 
(17) and (18), and store them. If they, respectively, are less 
than the Cosine and Sine of -48 degrees, equate them, 
respectively, to those of -48 degrees, 

d. Form (m+k) x (m+k) size matrices [Y8] and [YV] of (1) 
and (2) respectively each in a compact storage exploiting 
sparsity. The matrices are formed using (22), (23), (11), 
and (24). In [YV] matrix, replace diagonal elements 
corresponding to P V-nodes by very targe value (say, 
10.0 10 ). Factorize [Y6] and [YV] using the same ordering 
of nodes regardless of node-types and store them using 
the same indexing and addressing information. 

e. Compute residues [AP] at PQ- and PV-nodes and [AQ] at 
only PQ-nodes. if all are less than the tolerance (e), 
proceed to step-n. Otherwise follow the next step. 

f. Compute the vector of modified residues [RPJ using (20) 
for PQ-nodes, and using (9) and (19) for PV-nodes. 

g. Solve (3) for [AO] and perform an update, [6] = [6] + [a6]. 

h. Set voltage magnitudes of PV-nodes equal to the specified 
values, and Increment the iteration count ITRP=iTRP+l 
and r=(ITRP+ITRQ)/2. 

i. Compute residues [AP] at PQ- and PV-nodes and [AQ] at 
PQ-nodes only. If all are less than the tolerance (r.), 
proceed to step-n. Otherwise follow the next step. 

j. Compute the vector of modified residues [RQ] using (21) 
for only PQ-nodes. 

k. Solve (5) for [AV] and update PQ-node voltage 
magnitudes using [V] - [V] + [AV]. While solving (5), 
skip all the rows and columns corresponding to PV- 
nodes. 

I Calculate reactive power generation at PV-nodcs and tap 
positions of tap changing transformers. If the maximum 
and minimum reactive power generation capability and 
transformer tap position limits are violated, implement the 
violated physical limits and adjust the loadflow solution- 
is Increment the iteration count ITRQ=ITRQ+1 and 

r~0TRP+ITRQ)/2, and Proceed to stcp-c. 
n. From calculated and known values of voltage magnitude 
and voltage angle at different power network nodes, and 
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tap position of tap changing transformers, calculate power 
flows through power network components, reactive power 
generation at PV-nodes, and real and reactive power 
generation at the slack-node. 

IV. TEST RESULTS 

The proposed method is tested with a number of well-/ill- 
conditioned sample systems. Tt is also tested for the 
convergence sensitivity to branch R/X ratios. Two series of 
tests are generated* one by lowering the base case branch 
reactance and ihc other by raising the base branch resistances. 
Both scaling of a parameter of each of the branches 
simultaneously (uniform scaling) and scaling of a parameter of 
only one branch at a time (non-uniform scaling) are used. The 
latter generates more realistic difficult cases of a wide R/X 
ratio disparity among branches incident to a node. Wide 
disparities are also created at the terminal nodes of branches of 
zero resistances with the uniform scaling. 

When an increased R/X ratio is created by lowering 
reactance and leaving resistance the same, the branch 
impedance decreases, and improves voltage regulation for the 
receiving terminal node of a branch. The situation can be 
considered to be dominaatly simulating R/X ratio test 
condition [6]. A similar R/X ratio can be created by raising the 
branch resistance and leaving reactance unchanged. This 
increases the branch impedance and causes poor voltage 
regulation for the receiving terminal node of a branch. For a 
five-node network under uniform R -scaling test condition of 
(7.lR+jX). the minimum voltage magnitude in the solution is 
found to be about 0.5 p.u, and the system is very close to its 
voltage stability limit Larger scaling factor caused divergence 
.of the solution process. Such test conditions can be said to 
have simulated system non-linearity, which dominantly 
influence the convergence. 

All the tests are performed in single precision on personal 
computer using slack-start initialization. Results of the tests 
arc given in terms of the number of iterations required to 
converge to 0.01 Mw/Mvar mismatch APp/AQp unless another 
speeilic tolerance is mentioned. 

Programming of the polar co-ordinate formulations of the 
loadtlow methods studied implements rectangular power 
mismatch calculations involving storage of voltage magnitude 
vector and evaluation of trigonometric functions, sine and 
cosine of AO. in the iteration process [9]. Results are provided 
for trigonometric functions evaluation in all iterations. 
However, studies on different systems have shown that it is 
sufficient to evaluate trigonornctrie functions only in the first 
iteration with A8 assumed zero in subsequent iterations. 

Since off diagonal elements are made of branch 
susceptanccs and admittance magnitudes in SSDL-YY 
method, calculation of gain matrices involve taking only a 
square-root for each branch except those connected between 
two PV-nodes and a PV-nodc and the slack-node, and each 
diagonal clement takes only 4-mu Implication, a division, and 
a subtraction operation for I s 1 at which value all the method 
of this paper gives the best possible convergence. Also h 
takes only two multiplications and an addition operation more 
than the GFDL method at a PQ-nodc in each iteration. This 
overhead of super decoupling is negligible and can be ignored. 
Therefore, only iteration counts are used as the basis of 



comparison of proposed SSDL-YY method with the SSDL- 
BX and CFDL methods. 

Presently the only reasonable comparison in this paper is 
that between the GFDL method, its advanced version 
manifested as SSDL-BX method, proposed SSDL-YY method 
and FSDL [6] method. Results of Tables- 1, Table-2, Table-3, 
Table-5 and Table-6 provide this comparison. Clearly the 
SSDL-YY method is to be recommended. The comparison can 
be summarized as follows: 

(a) Relatively better performance of all the methods presented 
in this paper could be achieved for values of J= 0, 1, 2. 
Other values of I caused increased iterations to obtain 
converged solution. The best performance could be 
obtained for 1=1 for all of the methods, and therefore, test 
results of SSDL-BX and the proposed SSDL-YY methods 
are provided only for I«1 , 

(b) Table- 1 gives the convergence results for GFDL, SSDL- 
BX, SSDL-YY, and FSDL methods under uniform R- 
seale and X-scalc test conditions. The comparisons in (c), 
(d), (e) and (0 are based on results of Table- 1 - 

(c) The SSDL-BX method takes maximum of 6-5, 2-1, 4-4, 
9-8, and l-t iterations less than the GFDL method for 5-, 
14-, 26-, 57-, and 118-node systems, respectively, under 
uniform R-scaling test conditions, whereas it Uikes 
maximum of only l-l and 0-1 Iterations more for 14-node 
and 1 18-nodc systems respectively. SSDL-BX converges 
better for all systems tested over complete range of 
uniform X-scaling test conditions. Moreover it can be 
seen from Table- 1 that SSDL-BX method converges in 6- 
instances when the GFDL method £ails to provide 
solution, whereas it fails only in two instances when 
GFDL method provides solution. Therefore, SSDL-BX 
method is more efficient and reliable than GFDL method. 
This fact is further established from the anarysis of the 
number of test case given in Tablc-3. 

(d) The factor as determined in (19) is restricted to 0 75 to 
begin with, and tt is changed to 0.6 when the average of 
all the less than 1 0 value at PV-nodes is less than 0.6 for 
all versions of SSDL method in this paper. This iwo level 
restriction to K p is also applied to FSDL method. II can be 
seen that thus modified FSDL method tested for this paper 
itself is the substantial improvement over originally 
presented by this author in [6] particularly for 57-nodc 
system and under uniform X-scale factors 0.1 and less. 
Now, from the results of Table- 1, it can seen that SSDL- 
YY method takes maximum of 1-1, l-I, 3-4, 1-2, and 1-1 
iterations less than FSDL method for 5-, 14-, 26-, 57- and 
118-node systems respectively under uniform ll-scaling 
test conditions. The only exception is the uniform R- 
scaling test condition (3.5R+jX) where SSDL-YY method 
takes 1-1 iteration more in case of 14-nodc system. 

(c) SSDL-YY method takes maximum of 1-1, 1-2, 3-4, l-l, 
0-1 or 1-0 iterations less than FSDL method for 5-, 14-, 
26-, 57-, and 118-node systems respectively under 
uniform X-scaling test conditions for scale factors greater 
than 0.2, whereas it takes only l-l iteration more in case 
of 57-node system. For scale factors less than 0.2, SSDL- 
YY method takes maximum of 2-2, .V4, and 7-7iterations 
less for 5-,26~, and 1 18 -node systems respectively, where - 
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Toole- 1 : Iterations far convergence to 0,01 Mw/Mvar (Uniform scale factors) 



Test 5-nudc System [ 1 0] 

Foctor GFDL SSDL FSDL 

______ BX YY 



14-node System f J I] 
GFDL SSDL F$DL 
BX YY 



26-nodc Sy5(cmll2) 
Cl-DL SSDL FSDL 
BX YY 



57-nodc System [II J 



1 1 g-nodc Sytem 



CFDL SSDL FSDL QFOL SSDL FSDL 



Incrcwwd 

R+jX 
L3R+JX 

2R+jX 
2.5R+JX 

3R+iX 
3.3R+jX 

4JR+jX 

dftijx 

MR >jX 
7RljX 
7JR+JX 



n On- linearity condition tests generated by application of unitbmi R^c;tlc factors 



4-3 
4-4 

4- 3 
$4 

5- 4 

6- 6 
8-7 
9.8 
0-3 

10-9 

n-io 

13-13 



3- 2 
3^3 

4- 3 

5- 4 

6- 5 
6-5 

6- 5 

7- 6 

8- 7 



16-15 10.10 
ne nc 



3-2 

3- 3 

4- 3 
4^ 
4-3 
4-4 
4-4 

4-3 

4- 4 

5- 4 

5- 5 

6- 6 
7^7 
13-13 
nc 



4-3 
4-3 
4-3 
4-4 
4-4 

4- 4 

5- 4 
5-4 
5-5 

5- 5 

6- 6 

6- 6 

7- 7 
13-13 

nc 



4-3 


3-3 


3-3 


4-3 


8-8 


4-4 


5-4 


8-8 


5-4 


4-3 


4-3 




8-8 


4-4 


5-4 


8-8 


6-5 


4-4 


4-4 


5-4 


M 


5-5 


6-5 


9-9 


6-5 


5-4 


5-5 


6-5 


8* 


6-5 


6-5 


9-9 


6-5 


6-5 


5-5 


6-6 


9-8 


6-5 


6-5 


9*9 


7-6 


7-6 


7-6 


7-6 


10-9 


7-6 


8-7 


9-8 


8-7 


9-8 


9-8 


8-7 


12-12 J2-I2 


14-13 


16-15 


13-12 


14.13 


13-12 


13-12 


nc 


ne 


nc 


nc 


nc 


nc 


nc 


nc 













BX 


YY 






BX 


YY 




6-6 


7-6 


7-6 


7-6 


4-4 


4-3 


4-3 


4-4 


7-7 


7-6 


7-6 


7-6 


5-4 


4-4 


4-3 


4-4 


9-9 


7-6 


7-6 


8-7 


5-5 


5-5 


4-3 


5-4 


12-11 


8-7 


8-7 


8-8 


6-6 


6-5 


5-5 


5-5 


14-14 


10-9 


9-8 


9-9 


7-6 


6-6 


6-5 


6-5 


18-17 12-11 


10-9 


11-11 


7-6 


7-7 


H-7 


8-7 


23-22 


14-14 


12-tl 


13-12 


Q-Q 


9-9 




10-9 


26-25 


nc 


18-17 


18-17 


15-15 


14-14 


14-14 


14-14 


ne 




nc 


rtC 


nc 


nc 


nc 


ne 



increased 


brunch K/X ratio teats generated by application of uniform X-scaIc factors 
















R4jX 


4-4 


3-3 


3-3 


4-3 


5-4 


4-3 


4-3 


44 


8-8 


4-4 


5-4 


8-8 


7-7 


7-6 


7-6 7-6 


5-4 


4-4 


4-3 4-4 


R+J.9X 


4-4 


4-3 


3-3 


4-3 


5-5 


4-3 


^ 


5-4 


8-8 


4-4 


5-5 


8-8 


7-7 


6-6 


6-5 7-6 


5-4 


4-4 


4-4 5-4 


R*j.8X 


4-4 


4-3 


3-2 


4-3 


6-5 


4-3 


4-4 


5-4 


R-S 




5-5 


8-8 


7-7 


6-5 


6-5 7-6 


64 


4-4 


5-4 5-4 


R4j*7X 


4-4 


4-3 


3-2 


4-3 


6-5 


4-3 


5-4 


5-5 


7-7 


5-4 


6-5 


8-8 


8-7 


6-5 


6-5 6-6 


5-5 


5-4 


5-4 5-4 


Rfj,6X 
Rrj5X 


3-4 


4^3 


3-2 


4-3 


6-5 


4-3 


5-5 


6-5 


7-7 


5-5 


6-5 


7-7 


8-8 


6-5 


7-6 6^5 


6-5 


5^ 


5-4 54 


5^1 


4-3 


3-2 


4-3 


7-6 


4-3 


6-6 


7-6 


7-7 


6-5 


6-S 


7-7 


9-8 


6-6 


7-7 7-6 


6-6 


5-5 


6-5 6-5 


R+j.4X 


5-3 


4-3 


3-2 


4-3 


8-7 


4-3 


6-6 


7-6 


7-6 


6-5 


6-5 


7-7 


9-8 


7-6 


7-7 7-7 


7-7 


5-5 


6-5 6-5 


R+jJX 


6-5 


4-3 


3-2 


4-3 


8-7 


5-4 


6-5 


6-5 


8-7 


6-6 


6-5 


7.6 


10-10 


9-8 


8-7 8-8 


8-8 


6-6 


R-R R-R 


K i j,2X 


7-6 


4-3 


4-3 


4-4 


9-8 


5-5 


6^5 


7-7 


9-9 


7-6 


7-6 


7-7 


13-12 


12.11 


10-9. 9-8 


10-10 6-6 


9-8 9-8 


lUjJX 


y.y 


4-3 


5-4 


6-5 


12-1 1 


8-7 


9-8 


8-8 


12-11 


8-7 


13-12 


15-15 


30-29 


nc 


14-13 14-13 


nt 


19-19 


17-17 24-24 


R*j,09X 


|0-y 


4-3 


5-5 


6-5 


12-M 


9-8 


9-8 


8-8 


13-12 


8-8 


15-14 


IS-1H 


nc 




18-17 16-15 




27-26 26-25 24-23 


Rij.OSX 


12-11 


4-3 


5-5 


6-5 


12-11 


94? 


[f>-0 


9-8 


20-19 


10-9 


16-16 


15-15 






24-23 22-21 




nc 


«C nc 


R*j-07X 


1 J -10 


4-3 


5-5 


7-6 


13-12 


10-10 


10-10 


0-9 


ne 


15-15 23-22 


22-21 






34-33 34-33 








R rj»06X 


11-10 


5-4 


6.5 


7-7 


21-21 


14-13 


12-11 


9-9 




nc 


IK 


nc 






nc nc 








RHj.OSX 


13-12 


5-4 


6-5 


7-7 


ne 


nc 


26-25 21-20 






















R+j,04X 


nc 


6-5 


6-5 


7-7 






nc 


nc 






















R<j\03X 




7-6 


7-6 


8-7 


































6^5 


7-7 


9-8 






























Rij.OlX 




no 


7-7 


9-9 































ne : Not converged in 1 00 iterations or diverged 



as it takes maximum of 5-5, 1-1* 2-2, and 2-2 iterations 
more for 14-. 26-, 57% and 1 18-nodcs systems, 

(0 From (d) and (c) in the above it can be said that SSDL- 
YY method is overall better convergent than modified 
FSDL method tested for this paper, which in itself is an 
improvement over the one originally presented in [6]. 

(g) Toble-2 gives average number of iterations. Many 
load flows are run for each non-uniform scale factor. 
Load/low cases are generated by applying the same factor 
one at d time to different branches of a system. The 
average for 5*node system arc over 7-cascs of 7-branches 
for each scale factor. The averages arc over 32-cases of 
32-branchcs for26-node system, and they are over 10- 
cascs of first I (Mines for U8-node system. For each of 
the systems studied, individual R/X ratios are raised by 
the highest factor of 20 for non-uniform R-scaling and 
lowest factor of 0,03 for non-uniform X-scaling. It can be 
seen that a reduction in the average iterations is achieved 
for the SSDL-YY method compared to both SSDL-BX 
and FSDL methods. Average iterations for GFDL method 
arc not given as it is established in (a) above that SSDL- 
BX method is better convergent and more reliable than 



GFDL method. The reduction can be as much as about 
230%, 124%. and $7% against SSOL-BX method under 
non-uniform R-scaling test conditions, and 320%, 48%, 
and 130% under non-uniform X-scaling test conditions 
for 5-, 26-, and 1 18-node systems respectively. Moreover, 
the reductions in average iterations for SSDL-YY method 
compared to FSDL method are about 1 2%, 80%, and 14% 
under non-uniform R-scaling test conditions, and about 
13%, 80%, and 14% under non-uniform X-scaling test 
conditions for 5-, 26-, and 1 1 8-nodc systems respectively. 
However, slight increase of about 3% and 5% in average 
iterations by SSDL-YY against FSDL method is Observed 
under non-uniform X-scaling test conditions for S-, and 
1 1 8- node systems, 
(h) Table- 3 gives the analysis of total number of cases used 
to produce results of Tablc-2 for different test power 
systems. For each branch, seven non-uniform R- and X- 
scale factors other than the base case of Table-2 are con- 
sidered. Thus the total number of cases are 49, 140, 224, 
287, 70 and 70 for 5-nodc system of 7-branchcs, 14-nodc 
system of 20-branchcs, 26-nodc system of 32- branches, 
30-nodc system of 4 1 -branches, first 10-branches in 57- 
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TilMc-2: Average iteri twins for convergence to OjOl Mw/JVfvar (Non-uniform scale factors) 



T«l conditions 



Number of system nodes 



NX 
IDR+jX 

R*jO IX 
A+jO 05X 
R*j0.03X 



SSDL-BX 

3.00-3.00 
529-471 
S.Hrt-a.29 
11.57-11.00 
1371-1357 



5 -node 
KSDL 

4.00-3.00 
4.29-3.K6 
4.71^1.43 
4.71-4.57 
5.00-437 



3.00-3.00 4.00-3.00 

9.14-R.29 5.00-4.14 

10.14-15,29 6,00-5.57 

26.14-25,29 6.43-5.86 



SSDL-YY 


SSDL-BX 


26-nodc 
FSDL 


SSDLrYY 


SSDL-BX 


118-nodc 
FSDL 


SSDL-YY 


3.00-3.00 


4.00^.00 


8.00-8.00 


5,00-4.00 


4,00-4.00 


4.00-4.00 


4.00-3.00 


4.00-3.43 


4,83-4.50 


8,19-8.1-0 


5.13-4J8 


5.00-4.30 


5.104.40 


4.9TM.50 


4.57-4.14 


6.56*6.3 1 


8.34-8.34 


5.44^.78 


6.80-6.30 


5.IO-4.70 


5.10-4.70 


4.43-4.00 


8.41-H.19 


8.56-8.53 


5.72^.94 


7.80-7.50 


530-4.80 


5.00-4.60 


4.43-4.14 


n. 78-1 1.63 


8.88-8.81 


5.91-5.19 


8.90-8.60 


5.10-4.90 


4.90-4.60 


3-00-3,00 


4,00-4.00 


8.00-8*00 


5.00-4.00 


4.00-4.00 


4.00-4.00 


4.00-300 


5.29-4.5? 


5.75-5,28 


8.38-8.38 


6.00-5.34 


8.40-7.60 


(1.70-6.00 


6.80-6.30 


6.00-557 


R.00.7J9 


8.63-8.44 


6.72-6.03 


13-50-12.60 


R.4On8.30 


8.00-7.40 


6.57-6.00 


10.81-10.31 


9.31-9.00 


731-6.63 


20.10-19.30 


9.90-9.40 


9.10-8.40 



node system, and first 10-branches in 118-node system 
respectively. The results clearly indicate that the convergence 
performance and the reliability of obtaining solution is the best 
for the SSDL-YY method among the methods tested for this 
paper, followed by FSDL, and SSDL-BX methods. This fact is 
also established from the results of Table-I. 



Tablc-3; Convergence results for cases of non-uniform scaling 

Mclhotl $ysicm Number of cases solved in N No. 

Number iterations where N is Of 

Node* N£5 $<N5lO KXN520 20<Ns50 N>50 D1V Coses 



GFDL 5 


5 


22 


14 


8 


0 


0 


49 


14 


44 


43 


38 


10 


2 


3 


140 


26 


0 


170 


45 


5 


1 


3 


224 


30 


63 


94 


95 


24 


2 


9 


287 


57 


0 


50 


14 


6 


0 


0 


70 


118 


11 


35 


17 


7 


0 


0 


70 


SSDL-BX 5 


5 


21 


15 


8 


0 


0 


49 


14 


44 


44- 


37 


10 


2 


3 


140 


26 


109 


68 


41 


4 


2 


0 


224 


JO 


71 


86 


95 


25 


2 


8 


287 


57 


0 


52 


13 


5 


0 


0 


70 


MS 


12 


35 


19 


4- 


0 


0 


70 


r*$DL 5 


29 


20 


0 


0 


0 


0 


49 


14 


105 


30 


1 


0 


1 


3 


140 


26 


0 


209 


15 


0 


0 


0 


224 


30 


212 


67 


3 


0 


0 


5 


287 


57 


0 


63 


4 


1 


0 


0 


70 


118 


28 


36 


6 


0 


0 


0 


70 


SSDL-YY 5 


30 


19 


0 


0 


0 


0 


49 


14 


106 


29 


1 


0 


I 

1 


l' 

J 


140 


26 


139 


77 


8 


0 


0 


0 


224 


30 


230 


49 


3 


0 


0 


5 


287 


57 


0 


67 


2 


1 


0 


0 


70 


118 


32 


32 


6 


0 


0 


0 


70 



(i) Tablc>4 gives number of iterations talcen by the SSDL- 
YY method for different power mismatch requirements 
under three uniform scaling test Conditions including the 
base-cose of (R+jX). The number of iterations rises 
smoothly in line with the tolerance. 

(j) Toblc-5 gives number of iterations taken by different 
moihod for 13-node and 43-node ill-conditioned systems, 
ft should be noted that better convergence by SSDLrYY 



Table-*: Iterations by SSDL-YY for different tolerance* 



Uniform Tolerunrc 
Scaling test in Mw/Mvttr 



Number of system nodes 



Condition 




5 


14 


30 


57 


118 


2Ri jX 


1.000 


2-1 


3-3 


3-2 




3-2 




0.100 


3-2 


4-4 


4-4 


6-5 


4-3 




0.010 


4-3 


5-5 


5-4 


8-7 


5-5 


R+jX 


1.000 


2-1 


2-2 


3-2 


4-3 


3-2 


o.ioo 


3-2 


3-2 


3-2 


5-4 


3-3 




0.010 


3-3 


4-3 


4-3 


7-6 


4-3 


R+J0.2X 


L00O 


2-1 


4-3 


4-3 


5-4 


5-4 




0.100 


3-2 


54 


5*5 


7-6 


7-6 




0.010 


4-3 


6-5 


8-7 


10-9 


9-9 



method could be achieved by tuning rotation angle <I> p to - 0° 
for both the systems. The better convergence by FSDL 
method could be achieved by tuning O p to - 0° for 1 3-node 
system and by tuning <D p to-90 c for 43-node system. 



Tftble-5: tftmlioosfor convergence to 0.01 Mw/Mvar 


System 
No. Of 
Nodes 


l_ Method 


SSDL-YY 


FSDL 


SSD1^B> 


GFDL 


13 [131 
43fl31 


-48° -0° -90° 


-36° -0° -90° 


6-5 
14-14 


4-4 
18-17 


9-9 5-4 24-24 
20-19 15-14 20-19 


6-5 5-4 20-19 
15-14 29-28 12-11 



(k) Table-6 gives number of iterations taken by different 
methods to converge for each of Ihe 2-small, and 6- radial 
distribution systems. It can be seen that SSDJUYY 
method always converge in less number of iterations than 
those by FSDL method. The SSDL-BX method also 
converges in less number of iterations than by GFDL 
method in all except 40- node system in which case it 
takes almost double the iterations. GFDL method 
diverges in case of 31 -node system indicated by double 
star (**). However, SSDL-YY method converges in 
comparable to substantially less iterations than SSDL-BX 
method. Therefore; again SSDL-YY method turns out to 
be the winner among all. 
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Tftblo-tf: llcruHnm for convergence to 0.01 Mw/Mvar 


MfiTHOO 


Number of System Nodes 


3 


if 
.[»«] 


13 
1 151 


15 
JIT] 


22 
f ,5 T 


3! 


40 

[16) 


85 


SSDL-YY 


5-4 


6-5 


4-4 


2-1 


5-4 


7-7 


7-7 


4-3 


FSDL 


7-7 


7-f» 


5-4 


3-2 


R-7 


9-8 


9-8 


4-4 


ssouox 




6-5 


4-3 


3*2 


5^4 


2$-25 


28-27 


5-4 


CPDL 


7-7 


7* 


4-4 


4-3 


7-6 


*• 


14-13 


7-6 



V. FURTHER RESEARCH DIRECTION 

There are two rurthcr research directions. One is to study 
convergence pattern in the iterative solution process of the 
proposed method, and develop acceleration factors that can 
produce solution in maximum of 2 to 3 iterations. The other is 
to use -uniform and nonuniform P-scale and Q-scale factors 
for real power and reactive power loads in addition to uniform 
and non-uniform R-scale and X-SCale factors in training 
Artificial Neural Network (ANN) using the proposed method 
lo achieve very fast and almost closed form of solution, 

VI. CONCLUSION 

A generalised SSDL-X7C toadflow method is presented 
leading to SSDL-BX and the proposed SSDL-YY methods. 
The proposed method is characterized by gain matrices 
defined independent of rotation angles, rotation angles 
restricted to maximum of -48*, modification of real power 
mismatches at PV-nodes to correct for model inaccuracy, an 
efficient slack-start as initial guess solution, and the presence 
of transformed specified reactive power in the diagonal 
elements of ihc gain matrix [YV] of the Q-V sub-problem. 

The proposed SSDL-YY method is the most efficient and 
reliable among all the decoupled loadflow methods. The 
convergence of the method is demonstrated to be better than 
that of the SSDL-BX method for all normal R/X ratio systems 
to those having high R/X ratio branches. For systems of large 
R/X ratio branches, the proposed method can be more than 
twice faster than the SSDL-BX method, which is more reliable 
and better convergent than the GFDL method. It is also found 
lo be better convergent than the FSDL method. 

Attractive features of the proposed method are that it is not 
very sensitive to the restriction applied to the 
rotation/transformation angle, and it provides an elegant 
formulation for the diagonal elements of the gain matrix [YV] 
that suggest a mechanism lor their numerical value 
manipulation particularly useful when diagonal dominance 
issue arise in the presence of a capacitive series branch or 
excessive capacitive compensation at a node. 

Existing decoupled loadfllow programs can be easily 
upgraded to the SSDL-YY as is apparent from the description 
of the method. It is believed that all applications of the 
decoupled loadflow methods would benefit from the proposed 
method, and it would potentially be the preferred general 
purpose loadflow method. 

VII. APPENDIX 

A* Active and reactive power mismatch equations divided 
by V* on both sides and partial derivatives 



APp/Vp' - P/TV, 1 - (O w +gp)/Vp<»> 

- ^(C^CosGp, + B w Sin9 ro )V t /V p r, - l > (32) 

- Z (OpqSindn - B„Cose pq )V < /V p "- , > (33) 

•^APp/Vptyae,, = Zi (-G w Sm9p,+ Bp, Cwd^VjV™ (34) 

^(AQp/VpVdep . I (G„ Cose„ + B w Sin e N )V ( /Vp (, -'> (35) 

-ftAPp/Vp'ySe,, = (G„ Sine*, - Bp, Cosep,,)V < /Vp"- , > (36) 
-^AQ/Vp'yae,, = -(Gp,, Cos8 w + Op, SinG^V/Vp" 11 (37) 
-^APp/Vp'yaVp- iPp^/Vp^"- a-2)(Gpp+g p )A^ p (, -' ) 

- 0-1)1 (G w Cose (XI +Bp q Smep„)V < /Vp l (38) 

-3(AQ,/VpV3V p = T<yVvV w >+ (I^XBpp+b^yv, 0 ^ 

- (I-l)Z (OnqSinOpq - BpXosQp^V/vV (39) 

-d(^? ( /v p 1 ydv q - (Gpq Cos e M + b w sin e^yVp 0 ' 1 > (40) 

-^AQp/Vf/ySVi, = (Gpq Sin 8 pq - 6 W Cos e^yVp 0 ^ (41 ) 

In the fixed Jacobian methods, the Jacobian is to be 
determined only once at the start with an initial guess. 
Considering the guess solution to be the stack-start (V^6,for 
all the nodes), the expressions for the matrix elements reduce 
to the following: 

^ap ApVaa^Bpp/v ^ -3(A<yv 0 'y3e P = -o^/v™ 

^AP^pVaVpKlP^W^Gpp^T^Jgp)]^^ 
^AO p ^ p I )/^p=nQpW-(B pr H:i-2)bp)]^ a - 0 
-3(APp/Vp , )/3e q = -Bp/V^^ -3(APpA^p 2 )/3 V q = Gp^ Vp (M > 

Where, T^-oo, -2,-1. 0, 1,2 + ©©. 

The elements of the transformed gain matrices [Y6] and 
[YV] determined from the generalized Jacobian elements given 
in the above remain the same for any value of *T* except 
diagonal elements of the gain matrix [YV] of the Q-V sub- 
problem. The parameter h p that appears in the diagonal 
clement of the gain matrix [YV] can be expressed by (12) for 
different value of the raised power *T* for ihc generalized 
model SSDL-X'X\ and it is expressed by (24) for the 
proposed SSDL-YY model. 

The equations (12) and (24) provide elegant formulation for 
the diagonal elements of the gain matrix [YV] that suggest a 
mechanism for their numerical value manipulation particularly 
useful when diagonal dominance issue arise m the presence of 
capacitive series branch or excessive capacitive compensation 
at a node. It is always desirable to obtain a solution even if it is 
computationally expensive against not to have one. 

B. Transformed Power Flow Equations V 
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Real and imaginary parts of transformed power flow 
equations from [4] and [6] are given as: 



fV " vAGw-'+Sp') + VpKGp,' Cos V-Bp,' Sin O^V, 



(42) 

(43) 

(44) 
(45) 
(46) 



Qp' - -V^pp'+V) * VpZXGp,' Sin6 w - B„' Cose„)V, 

where, G M ' + jB^ * Y N exp (j(Q> n + 1 80* + Op)) 

Cpp* + j*V - Ypp exp Gl^pp + *p» 

gp* + jV - y P e^p G(<V + Op)) 
C Transformation of Branch Admittance 

The gain matrices in their transformed version arc 
asymmetrical because mostly different rotations are required 
to be applied at the terminal nodes of a bmnch [4}. However, 
symmetrical gain matrices can be obtained by the Haley and 
Ayrcs technique [5] of applying an average of rotations at the 
terminal nodes of a branch to the branch admittance, and it 
involves calculation of an inverse of the Tangent at each node 
ro know the rotation angle and the Sine as well as Cosine 
trigonometric functions of the average of the rotations ai the 
terminal nodes of a branch. 

The constant jrymmeirical gain matrices of the Decoupled 
Loadflow Models SSDL-YY, SSDL-XX, and SSDL-BX are 
independent of rotation angles. However, gain matrices of the " 
$SDL-X*X* and other models described in [18] are dependent 
on rotation angles as their elements are determined by 
transforming a branch admittance as shown in following steps 
in order to keep the gain matrices symmetricaL 

I > Compute: flyarctan (Gpp/Bpp) &O q ^ arctan (G^/B^ (47) 

2- Compute the average of rotations at the terminal nodes (p 
and q) of a branch; «P», = (<X> P + OJ/2 (48) 

Compare <P„ with the Limiting Rotation Angle (LRA) 
and let 0 W to be the smaller of the two: 

d>av ■ minimum (O avi LRA) (49) 

Compute transformed pq-th element of the admittance 

matrix; G OT *+jB w >ea (Cos O av + jSin o f¥ )(G w +jB Iiq ) 

Note (hat the transformed branch reactance is: 



3. 



4. 



(50) 



5. 



" B w 7 (<W 2 + Similarly, 
X^^BppVCGpp^H-Bpp' 2 ) 



(51) 
(52) 

In the description above X„' is the transformed branch 
reactance defined by equation (51) and B N * is the 
corresponding transformed element of the susceptance matrix. 
Gpp' and Bp/ are diagonal elements as given in (52). 

D* Loadflow Models organized for Global Corrections 

Tlw loadflow models described in this paper and all others 
can be organized to produce corrections to the initial estimate 
solution or global corrections in all iterations. This involves 
storage of the vectors of modified residues and replacing the 
equations (7), ft), (9), (50). (4) and (6) respectively by (53), 
(54). (55). (56), (57) and (58). Superscript '0* in equations 
(57) and (58) indicates the initial solution estimate. 



RPp r - E(AP P ') r + (GppVBpp'XAQp*) y(V/) ' + RP p < *'> 
RQn r « KA<V) r + (Gp P '/Bpp\)(AP p ') r yVp r + RQ (M) 



RP p r = AP p r /(KpVp r ) + RP p tM) 



RJY 



AP'/V/ + RP< M) 



ep r =ep o + A0p r 

V p r -V p ° + AVp r 



I. 



(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
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